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SUMMARY 


Two approximation concepts for combined thermal/structural design are 
evaluated. The first concept is an approximate thermal analysis based on the 
first derivatives of structural temperatures with respect to design variables. 
Two commonly used first-order Taylor series expansions are examined. The 
direct and reciprocal expansions are shown to be special members of a general 
family of approximations, and it is shown that for some conditions other mem- 
bers of that family of approximations are more accurate. Several examples are 
used to compare the accuracy of the different expansions. 

The second approximation concept is the use of critical time points for 
combined thermal and stress analyses of structures with transient loading con- 
ditions. It is shown that significant time savings may be realized by identi- 
fying critical time points and performing the stress analysis for those points 
only. This approach is used to design an insulated panel which is exposed to 
transient heating conditions. 


INTRODUCTION 

One of the major obstacles to the widespread use of automated structural 
design (member sizing) procedures is the need for large computational resources 
to perform repeated analyses of a structure throughout the design process. To 
alleviate this problem, it is now common practice during major portions of the 
design process to use approximations for response functions such as displace- 
ments, stresses, buckling loads, or flutter speeds. (See refs. 1 to 6.) Such 
approximations are often based on the first derivatives of the response func- 
tions with respect to the design variables. One obvious technique in the use 
of these approximations is to employ a first-order Taylor series (i.e., a lin- 
ear approximation) , herein called the direct expansion. An alternate approach 
is to approximate the desired function by a first-order Taylor series in the 
reciprocals of the design variables. Such a technique, herein called the 
reciprocal expansion, is exact for stresses and displacements of statically 
determinate trusses and has been shown to be more accurate than the direct 
expansion for other structures (ref. 5). 

Another useful approximation concept is that of constraint deletion 
(e.g., ref. 1). Constraints which are far from being critical are ignored 
during parts of this design process. Periodic updates are made of the number 
of constraints that are retained. For constraints which are functions of a 
parameter such as time (parametric constraints) , the concept of constraint 
deletion can be extended (ref. 7). When the dependence on the parameter 
is sufficiently smooth, critical values of that parameter (critical points) 
may be identified. The constraints for other values of the parameter may be 
ignored; however, the values of the critical points are updated periodically 
as the design changes. 



This paper is concerned with application of the two concepts, response 
approximation and critical points, to combined thermal/structural design. Both 
direct and reciprocal expansions of the response functions are investigated for 
accuracy in predicting temperature fields for steady-state and transient thermal 
problems. The time savings associated with the use of these expansions are also 
evaluated with the aid of two examples; one of a steady- state response and the 
other of transient thermal response. In the second example, use is also made of 
the second approximaticai concept (i.e., critical time points), and an investi- 
gation is made of the resulting accuracy and time savings. 

SYMBOLS 

Values are given in both SI and u.S. Customary Units. The measurements 
and calculations were made in U.S. Customary Units. 

a-| ,a2fa3 bar element areas, mm^ (in^) 

E Young's modulus. Pa (psi) 

f (V) response function 

fotV) direct aH>roximation to function f (V) 
fR(V) reciprocal approximation to function f (V) 
g constraint function 

g^ critical constraint function 

k thermal conductivity, W/m-^ (Btu-in/sec-f t^-Op) 

m(V) mass of structure, kg (lb) 

n^ number of constraints 

n^ number of design variables 

n-f number of time points 

Njf/Ny,NjQ, applied loads, kN/m (Ibf/in.) 

Q heat rate, W/m^ (Btu/in^-sec) 

t thickness, mm (in.) 

t-|,t2,t3,t4 plate element thicknesses, mm (in.) 

tQ,t45,t9o lamina thicknesses, mm (in.) 
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T temperature, (^P) 

Tg allowable temperature, {^) 

T^q transient, surface equilibrium temperature, (^) 

Vi ith design variable 

Vj^i offset value of ith design variable 

V vector of design variables 

a coefficient of thermal expansion, per (per 

p mass density, kg/m^ (Ib/in^) 

T time, sec 

Tc critical time, sec 

Subscripts: 

max maximum 

o initial value 


APPROXIMATION CONCEPTS 

First-Order Taylor Series Expansions of Response Functions 

The design problem considered herein is to find a vector of design vari- 
ables V (with components Vj; j - 1, 2, • • ., n^) that minimizes the mass 
m(V) of a structure and/or its thermal protection system subject to a set of 
constraints 


gi(V,T) > 0 


(i = l, 2, • • •, n^y 0 < t < (1 ) 


where x is time. The constraint functions g^ represent design limits on 
response functions such as displacements, stresses, or temperatures. Many 
algorithms for resizing the structure require evaluation of the derivatives 
of the response functions with respect to the design variables. These deriv- 
atives may be used to obtain approximations for the response functions which 
can replace the full analysis during parts of the design process. One approxi- 
mation for a response function f(V), which is based on the value of the func- 
tion and its derivatives at a design point Vor is a first-order Taylor series 
expans ion 
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af (Vo) 


(2) 


fu(V) = f(Vo) + 


j=i 


3vh 


(Vj - Voj) 


where Voj are the components of Vq. This approximation is herein called 
the direct expansion. An alternate expression (e.g.r refs. 1 to 3, 5f and 6), 
herein called the reciprocal expansion, is a first-order Taylor series in the 
reciprocals of the design variables. It can be written as 


V 3f(Vo) Voj 

fR(V) = f(Vo) + > (Vj 

4^ 3V-; v^ 


- Voj) 


(3) 


The reciprocal expansion was originally proposed for approximating stresses 
and displacements in structural analyses. For statically determinate trusses, 
stresses and displacements are indeed linear functions of the reciprocals of the 
cross-sectional areas of the members. For statically indeterminate structures, 
the reciprocal expansion is not exact, but studies (ref. 5) have shown it to be 
superior to the direct expansion for a number of problems. The reciprocal 
expansion also has the added attraction that it is usually more conservative 
than the direct expansion in estimating response functions such as stresses, 
displacements, or temperatures (ref. 2). If f(V) is the temperature at a 
point and V is a vector of insulation thicknesses, it may be expected that 
3f(V)/9vj < 0. Since all design variables are positive, it is easy to check 
from equations (2) and (3) that fR(V) > f£)(V) . 

The reciprocal expansion predicts an infinite value for f (V) when any 
one of the variables Vj goes to zero. It is, therefore, especially suitable 
for approximating functions, such as stresses and displacements in a stati- 
cally determinate truss, which go to infinity when the member areas go to zero. 
Approximation of stresses in redundant structures (or temperatures in situations 
where multiple heat paths exist) is slightly different in that several variables 
must vanish simultaneously for f (V) to become infinite. 

Some understanding of the relative performance of the direct and reciprocal 
expansions for a redundant structure may be gained by an examination of the 
following simple example: 


n 

f(V) = 

1 0v-| + V 2 


(4) 


This example illustrates a redundant situation where one of the design variables 
has a much larger effect than the other. Assume that f(V) is approximated 
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by using its value and derivatives at Vq] = Vq 2 1.0. The direct and recip- 
rocal expansions for several other values of v-j and V2 contained in 

table 1 • Two conclusions may be drawn from table 1 • First, the reciprocal 
expansion is better for approximating f(V) for changes in the dominant design 
variable V] but not as good for approximations with changes in the other design 
variable V2« Second, the reciprocal expansion is much better with changes in 
the design variables which reinforce each other (v] = V 2 = 2) than with changes 
that cancel each other (vj = 0.5, V 2 = 6) . 


Generalized Approximation 

The fact that the function f(V) becomes very large only when several 
variables go to zero simultaneously can be accounted for by the following more 
general form of the approximations 


n,; 


V ^oi " Vjni 3f 

f (V) s f(Vo) + > (Vi - Voi) — (Vo) 

Vi - vmi 3vi 

1 = 1 


(5) 


where v^^i is the ith component of a vector of constants chosen to account for 
the fact that as one variable goes to zero the function remains finite because 
other variables may have nonzero values. The direct and reciprocal expansions 
are special conditions of equation (5) with the direct expansion correspond- 
ing to v^i -oo and the reciprocal expansion corresponding to v^^i = 0. 
Because the approximation becomes unbounded when vi = Vj^ji, Vj^i may be chosen 
as the value of Vi which drives the exact function f (^0 to infinity if all 
other variables Vj[ remain constant (i.e., Vj = Voj where j ^ i) . 

For the function of equation (4) with Vqi = Vq 2 = ^ / such a consideration 
leads to the choice of Vj^i = -0.1 and v^2 " -10. The last column in table 1 
shows that the general approximation for this choice is generally superior to 
the other two, in the range considered. In general, it may not be possible to 
pick optimum values of v^^i, but small negative values of may be used to 

avoid large errors in the reciprocal expansion when one variable goes to zero. 


Critical Time Points 

The use of response approximations applies to designs under either steady- 
state or transient thermal constraints; however, a second approximation concept 
may also be useful in transient problems. For transient thermal loads, the 
thermal/structural design problem is characterized by constraint functions (see 
eq. (1)) which must be enforced for an entire time interval. Numerically, this 
means that a set of closely spaced time points Xj, where j=l,2, ..., n^, 
must be chosen such that serious constraint violations at intermediate time 
points are unlikely. Equation (1) may then be replaced by 
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= 2f • • \ 

( 6 ) 

j — If 2f • • »f 


In most practical applications, the product n^n^ will be a very large number 
and carrying out the design process with this full complement of constraints 
may become prohibitively expensive. 

Instead of monitoring a constraint at a set of n^ time points, it may 
be monitored only at its most critical points (ref. 7). This concept is 
explained by figure 1 which schematically shows the variation of a constraint 
function with time. The times of the local minima of the constraint function 
(A, B, and C in fig. 1) are herein called critical time points. Instead of 
monitoring the constraint function at all times, it is monitored only at those 
times. However, the critical time points drift as the design changes, and the 
constraint must be calculated at many time points to accurately locate the 
critical time points. The usefulness of the critical-time-points concept 
derives from the fact that, for portions of the design process, the drift of 
the critical time points may be neglected. For these portions of the design 
process, the constraint function has to be evaluated for only a relatively 
small number of time points. 

To prove that the drift in the location of the critical points may be 
neglected to a first-order approximation, it is assumed at first that a con- 
straint function g(V,T) has a single local minimum in the time region of 
interest. Equation (1) for the constraint gi(V,x) is replaced by 

9ci(V) = 9i(V,Tc) = min{gi(V,T)} ^ 0 (7) 


where Xc is the time for which g^ is most critical (the critical time 
point). The partial derivative of g^^ with respect to a design variable Vj 
is 


8gci(V) 3gi(V,Tc) 9gi(V,Tc) 3 tc 

= + ( 8 ) 

3vj 3vj 3 t 3vj 


For interior critical points (A and B in fig. 1), 3gi/3T = 0, and for boundary 
minimum points (point C) , 3Tc/3vj =0 if 3gi/3x 0. Therefore, 


3gci(V) 3gi(V,Xc) 

3vj 3vj 


(9) 
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Equation (9) shows that the drift in the critical time point may be neglected 
to a first-order approximation. Drift may be neglected not because it is small 
but because it does not affect the first derivatives of the constraint 9ci* 1^^ 
practice, the critical time points are calculated periodically and are frozen 
between updatings. However, the concept remains useful, even if the critical 
points require frequent updating, because of the savings in the derivative 
calculation. 

When the constraint function has more than one local minimum as in fig- 
ure 1 , each critical point must be assigned a separate constraint gci# other- 
wise, equation (7) may yield a discontinuous derivative of when the global 

minimum is switched from one local minimum to another. 


EVALUATION OF APPROXIMATION TECHNIQUES 

To judge the usefulness of an approximation concept, the question of what 
magnitude of change in design parameters results in unacceptable errors must 
be answered. Although it is desirable to obtain a general answer to this ques- 
tion, this did not appear possible in the present studies. Instead, the present 
work uses examples of both steady-state and transient thermal/structural design 
results to assess the usefulness of the two approximation concepts discussed in 
the preceding section. The steady-state example is a titanium panel with inte- 
gral aluminum bars, and the transient one is an insulated panel. 


Approximate Thermal Analysis of a Titanium Panel 

The accuracy of both the direct and reciprocal expansions for steady-state 
thermal analysis was investigated in reference 8 for a configuration which con- 
sisted of a titanium panel with aluminum bars as shown in figure 2. This con- 
figuration is representative of a class of structures in which one material 
satisfies strength requirements and the other acts as an efficient conductor 
to transfer incident heat to a heat sink. The sink is represented by the panel 
edge maintained at T = -18^ C (0^ F) . The incident heat load, material proper- 
ties, and remaining boundary conditions are also shown in figure 2. The accu- 
racy of the direct and reciprocal expansions for predicting the maximum temper- 
ature is shown in figure 3 for changes in the thickness of one-quarter of the 
plate and in the area of one aluminum bar (both variables are changed propor- 
tionately) , Initially the eight triangular plate elements were 12.7 mm 
(0,5 in.) thick and the four bar elements had areas of 645.2 mm^ (i.o in^) . 

The results shown in figure 3 indicate that the reciprocal expansion is 
more accurate than the direct expansion. The portions of the plate and the 
bars that were resized constitute the major path of heat conduction to the heat 
sink. When their thickness and area are reduced to zero, the temperature 
becomes unbounded so that the reciprocal expansion may be expected to be better 
than the direct one. The bias toward the reciprocal expansion can be carried 
to an extreme by uniformly increasing or decreasing the thickness and areas of 
all elements. When this is done, the error in the reciprocal expansion becomes 
minute, A different situation occurs when some of the design variables are 
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increased and others are decreased (a common occurrence in optimization prob- 
lems) . Performance of the expansions for such mixed changes in design variables 
is shown in figure 4. The percentage change in maximum plate temperature is 
shown as a function of the percentage of the total change in design variables 
fran the original design to the final design. This is indicated in the figure 
inset. For these conditions, the direct expansion is better. The reciprocal 
expansion does not perform as well because of the redundancy of the heat paths. 
This redundancy creates a situatiai similar to that of the simple example of 
equation (4 ) previously noted . 


Approximate Thermal Analysis of an Insulated Panel 

The quality of the approximate thermal analysis was investigated for tran- 
sient heating of the insulated panel shown in figure 5. A transient surface 
equilibrium temperature typical of a reentry heating trajectory, is 

applied to the outer surface of the insulation layer. The insulation protects 
a balanced, symmetric, graphite/epoxy composite panel with 0^, ±45^, and 90^ 
plies. The material properties of the composite panel and insulation are given 
in table 2. The thermal and stress analyses of the panel were carried out using 
a modified version of the computer program employed to obtain the results of 
reference 8. 

The dependence of the laminate temperature on insulation thickness is shown 
in figure 6 for one time point. Also shown are the direct and reciprocal expan- 
sions for that temperature based on derivatives calculated for a 1.52-mm 
(0.060-in.) laminate thickness and 102-inm (4.0-in.) insulation thickness. The 
reciprocal expansion is reasonably good over the entire range of insulation 
thicknesses, which is equal to the base value ±50 percent. The direct expansion 
is considerably less accurate. The better performance of the reciprocal expan- 
sion is understandable because the temperature becomes very high when the insu- 
lation thickness goes to zero. 

The dependence of the laminate temperature on laminate thickness is shown 
in figure 7. The temperature is much less sensitive to the laminate thickness 
than to the insulation thickness so that the errors for both approximations are 
smaller. For this set of calculations however, the direct expansion is better 
than the reciprocal one. The reason is that the temperature is not expected 
to become very high when the laminate thickness goes to zero because of the 
presence of the insulation (the same phenomenon was noted in the example of 
eq. (4) for changes in V 2 ) . It is also shown in figure 7 that an even better 
approximation can be obtained with Vj^ = -4.06 mm (-0.T6 in.). 

For the calculations illustrated in figure 7, the design variable was 
assigned to the laminate thickness. Often, design variables are assigned to 
the thickness of each ply. When one ply thickness is changed while others are 
kept constant, the reciprocal expansion will be even less accurate. This 
inaccuracy was demonstrated by repeating the calculations presented in fig- 
ure 7 with individual design variables assigned to the 0^, ±45^, and 90^ plies. 
The initial laminate thickness, 1.52 mm (0.060 in.), consisted of 0.5i mm 
(0.020 in.) for each ply direction. Then, the 90^- and ±45^-ply thicknesses 
were kept constant and only the 0^-ply thickness was increased. The reciprocal 
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expansion for these conditions is shown in figure 7 by the circular symbols. 

The direct expansion is not affected by assigning design variables in this 
manner . 

The curve with the circular symbols in figure 7 shows the danger of 
choosing Vj^^ = 0 (reciprocal expansion) for a highly redundant system such 
as a composite laminate. The choice of Vj^i is based on estimations of the 
ccanbined effect of all other design variables when v^ goes to zero. Based 
on the insulated panel example, the following recommendations can be made: 

1 . If the response function is not expected to become very large when v^ 
goes to zero, the direct expansion is appropriate and v^i should be chosen to 
be a large negative number . 

2. If the response function may be expected to become large but finite 
when V£ = 0, choose Vj^i to be a small negative fraction of a typical value 
of Vi* This is a safety measure to guard against excessive errors in the 
approximation when Vi changes by an order of magnitude. 


Use of Critical Time Points 

The approximation concept of following only critical time points was also 
checked for the insulated panel. Shown in figure 8 are typical changes in 
temperature histories as the design changes. Indicated in the figure is an 
appreciable drift in the critical time for the maximum temperature constraint. 
The same phenomenon is illustrated in figure 9 by showing the difference between 
following the maximum temperature, that takes into account the drift of the 
critical time point, and following the temperature at the time point which was 
maximum for the base design (insulation thickness = 101.6 mm (4.0 in.), laminate 
thickness = 1.52 mm (0.06 in.)). The two curves in figure 9 have the same 
derivatives at the base design point as predicted by equation (9). However, 
the difference between the two curves becomes substantial due to the drift of 
the time point of maximum temperature. There are indications in figures 8 and 
9 that the critical time points require updating more often than the approximate 
thermal analyses. (Compare figs. 6, 8, and 9.) 


APPLICATION TO OPTIMIZATION 
Optimization of a Titanium Plate 

The titanium plate with aluminum bars shown in figure 2 was designed for 
minimum mass under a constraint on the maximum temperature. A resizing 
algorithm based on optimality criteria (ref. 9) was used for optimization. 
Derivatives of the temperatures with respect to design variables were calcu- 
lated analytically. 

Mass iteration histories for this design problem are shown in figure lO. 
During the first few iterations, all design variables were reduced until con- 
straint violations occurred; thus, a reciprocal expansion may be expected to be 
more accurate than a direct expansion. Succeeding design iterations tend to 
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redistribute material among the structural elements to satisfy the design con- 
straints; therefore, a direct expansion may be expected to be more accurate 
for that phase of the optimization process. Accordingly, three separate cal- 
culations were made using approximate thermal analysis: in the first, a 

reciprocal expansion was used in both phases of the optimization process; in 
the second, a direct expansion was used in both phases; and in the third, a 
reciprocal expansion was used in the first phase and a direct expansion was 
used in the second phase. For the initial design, the plate thickness was set 
at 76 imi (3.0 in.) and the aluminum bars had areas of 9290 mm2 ( 14.4 in2) 
each. Derivatives were calculated for the initial design, and the temperature 
was approximated up to iteration 10. At that point, a new full analysis was 
used to generate a new approximation which was used for the rest of the design 
process. The decision to generate a new approximation at iteration 10 was 
based an full analysis results where changes in the design variables became 
mixed rather than uniform. Although the mass changed very little after about 
the 40th iteration, the optimality criteria required an additional 40 itera- 
tions to accurately satisfy the temperature constraints. 

The final designs and maximum temperature errors (which were obtained 
from a full analysis of that final design) are given in table 3. It can be 
seen that the direct expansion lacks sufficient accuracy to guide the opti- 
mization process completely to a converged solution. The reciprocal expan- 
sion also begins to diverge after about 30 iterations. This results in a 
design substantially different from the full analysis with 20 percent greater 
mass and a maximum temperature error of 49.9° C (89.8° F) . However, the use 
of the reciprocal expansion in the first phase of the optimization and the 
direct expansion in the second phase yields a final design very close to the 
full analysis design with only a i .1 percent mass difference and a maximum 
temperature error of 14.9° C (26.8° F) . Additionally, the approximate thermal 
analysis reduces the total optimization time by about 40 percent. The opti- 
mizaticxi algorithm was modified to automatically select the expansion concept 
based on whether the changes in design variables all have the same sign or 
different signs. The results obtained by such a modified algorithm were essen- 
tially the same as before. 


Optimization of an Insulated Panel 

The insulated panel shown in figure 5 was designed for minimum mass under 
temperature and stress constraints. The stress constraint was based on the 
Tsai-Wu failure criterion (ref. 1 O). The applied loads and the material 
strength allowables are given in table 2 . 

The optimizer for the programs for analysis and resizing of structures 
(PARS), described in reference 11 , was used for optimization. PARS is a 
collection of programs developed for the minimum mass design of structures 
subject to stress, displacement, and flutter constraints. The PARS optimizer 
employs the sequential unconstrained minimization technique (SUMT) , described 
in reference 12 , with a quadratic extended interior penalty function, described 
in reference 13. Newton's method with approximate second derivatives (ref. 14) 
is used to generate directions for one-dimensional searches for each uncon- 
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strained minimization. The program approximates the constraint functions dur- 
ing one-dimensional searches using values of specified by the user. 

In resizing the panel, new temperatures were approximated based on the 
values and derivatives of the initial design temperatures without any updating. 
That is, a full thermal analysis was performed only twice; once for the initial 
design and once for the final design. Exact calculations of the stresses from 
the approximate temperatures were performed once for each one-dimensional 
search. During each search, the stresses were approximated using the same 
values of used for approximating the temperatures. 

Starting with two initial designs, studies were carried out to test the 
accuracy and time savings associated with the two approximation concepts. The 
first initial design had a 101,6-mm (4,0-in,) insulation thickness and ply thick- 
ness of tg = tgo = 0,76 ram (0.03 in,) and t 45 = 1,52 mm (0,06 in,, combined 
thickness of +45^ and -45^ plies) , This is fairly close to the optimum insula- 
tion thickness, but far from the optimum laminate thickness. The second ini- 
tial design had a 203-mm (8.0-in.) insulation thickness and ply thickness of 
^0 “ ^45 “ ^90 “ 0.508 ram (0,02 in,), which is close to the optimum laminate 
thickness but far from the optimum insulation thickness. 

Performance using resp onse approximations The results of the optimiza- 
tion starting with the first initial design are summarized in the first three 
rows of table 4. For this design, both approximations compare favorably with 
the full analysis. This may be expected, because the final design differs 
from the initial design primarily in the laminate thickness which has a weak 
effect on the temperature. (See fig. 7.) The final design thicknesses deter- 
mined using the reciprocal expansion are substantially different from those 
obtained frcxn the full analysis. However, the mass is only 3.6 percent greater 
than that from the full analysis design. The maximum temperature error for the 
reciprocal expansicn final design is only 3.67^ C (6.6^ F) , which is smaller 
than the error for the direct expansion. This does not indicate that the recip- 
rocal expansicHi is superior to the direct expansion; rather, it indicates that 
the difference between the initial and final design thicknesses from the recip- 
rocal expansicai is less than that from the direct expansion- 

considerable time savings are achieved by using the approximate thermal 
analysis. The approximate thermal analysis requires less than 5 percent of the 
time required for a full analysis. However, the total optimization time can 
cnly be reduced to 30 percent of the time required for full analyses because 
of the time needed for the initial full analysis and for the operation of the 
PARS optimizer. 

The first three rows of table 5 summarize the results obtained with the 
second initial design. This time there was a large change in insulation thick- 
ness, and as a result, both approximations performed rather poorly. The direct 
expansicai yielded an infeasible design with a maximum temperature error of 
292^ C (525^ F) , The reciprocal expansion produced a much better design which 
is very close to the true optimum with a maximum temperature difference of 
66.7^ C (120^ F) . Based on the temperature error, the final design is, in 
fact, much closer to the true optimum than may be expected. This is because 
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the error peaks for noncritical time points. For the critical time points, the 
temperature error is an order of magnitude smaller. 

The results in tables 4 and 5 for the reciprocal expansion were obtained 
with Vmi = 0 for the ply thickness design variables. When the same calcu- 
lations were repeated for Vmi = -1.02 mm (-0.04 in.), there was little change. 
However, the changes in ply thickness from the initial design to the final 
design were moderate. When ply thicknesses are changed drastically during the 
design process, it becomes necessary to compensate for the thickness of other 
plies via Vjj,i. An example is an initial design with insulation thickness of 
122 mm (4.8 in.) and ply thickness of tQ = tgQ = 0.127 mm (0.005 in.) and ^ 

t 45 = 1.27 ram (0.05 in.). With Vmi = -1.02 mm (-0.04 in.), the final mass 
per unit area was 20.67 kg/m^ (4.233 Ib/f t^) , or 3.6 percent heavier than the 
true optimum, with a maximum temperature error of 7.4° C (13.3° P) . With ' 

Vjni = 0 mm, the final mass per unit area was 21.35 kg/m^ (4.372 Ib/ft^) , or 
7.1 percent heavier than the true optimum, with a maximum temperature error of 
16.4° C (29.6° F) . 

Based on these studies, it appears that it may be feasible to employ an 
approximate thermal analysis using the reciprocal expansion without updating 
the approximation. If the temperature error for the final design is large 
enough to put the optimality of that design into doubt, the optimization may 
be restarted from that final design. This was done for the final design given 
in table 5 for the reciprocal expansion. The restarted optimization required 
14.0 seconds of central processing unit (CPU) time to converge to the exact 
optimum in a few iterations. 

Performance using critical times .- The performance of the critical-time 
approach was evaluated next. Because of the drift in the critical time point 
discussed in "Use of Critical Time Points" and shown in figures 8 and 9, crit- 
ical times were updated at the beginning of each one-dimensional search. For 
each one-dimensional search, stresses were calculated only at two time points: 
the time of maximum stress ratio and that of maximum temperature. The results 
of the optimization using this concept are summarized in the last two rows of 
tables 4 and 5. It can be seen that additional savings of about 50 percent in 
computation time are realized without any substantial change in the details of 
the final design. 


CONCLUDING REMARKS 

Two approximaticxi concepts were evaluated for use in the design of thermal/ 
structural systems for minimum mass. The first concept was the use of response 
approximations, and the second concept was the use of critical time points. The 
evaluations were carried out by the use of steady-state and transient example 
problems. The following conclusions are therefore valid only to the extent that 
they apply to those examples. 

The first concept which was evaluated made use of approximate thermal anal- 
ysis. Two commonly used approximations, the direct and reciprocal expansions, 
were shown to be special members of a more general family of approximations 
characterized by a set of parameters v^i which controls the point where the 
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approximations become unbounded. The following conclusions were reached regard- 
ing the performance of the approximations. 

1. The reciprocal expansion is more accurate than the direct expansion when 
the temperature becomes very high as a single design variable goes to zero. 

2. If the temperature does not become very high when a single design vari- 
able goes to zero but as several design variables go to zero simultaneously 
(redundant system) , the reciprocal expansion provides a better approximation 
when these design variables either increase or decrease together. The direct 
expansion is better when some design variables are increased as others decrease. 

3. The reciprocal expansion is accurate enough for approximating temper- 
atures without needing updating by a full analysis for changes of up to 50 per- 
cent in the values of design variables during the design process. 

4. For a redundant system such as a composite laminate, the generalized 
approximation with a value of the parameter which compensates for the redundancy 
is found to be better than either the direct or reciprocal expansions. 

The second concept that was evaluated, the use of critical time points, was 
found to work well if the critical time points were updated periodically. For 
transient example problems, the combined use of the two approximation concepts 
resulted in an order of magnitude reduction in computation time. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
April 10, 1979 
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TABLE 1.- APPROXIMATIONS FOR f(V) = 11/(1 Ovi + V 2 ) BASED ON 
VALUES OF FUNCTION AND DERIVATIVES AT VqI = Vq2 = 1 


Vl 

V2 

f (V) 

Direct 

appr oxima t ion 
(eq. (2)) 

Reciprocal 
approximation 
(eq. (3)) 

Generalized 
appr ox ima t ion 
= -0*1^ ^m2 = -10 

(eq. (5)) 

1 .0 

1.0 

1 .000 

1.000 

1.000 

1 .000 

.5 

1 .0 

1.833 

1 .455 

1 .909 

1 .833 

2.0 

1 .0 

.524 

.091 

.545 

.524 

1 .0 

.2 

1 .078 

1 .073 

1 .364 

1 .078 

1 .0 

5.0 

.733 

.636 

.927 

.733 

2.0 

2.0 

.500 

.000 

.500 

' .441 

.5 

6.0 

0 

0 

. 

• 

0 

0 

0 

1.833 

1.521 
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TABLE 2.- MATERIAL PROPERTIES AND LOADING FOR INSULATED PANEL 


(a) Mechanical properties for graphite epoxy 


Property 

Symbol 

Units 

Value at - 

Room temperature 

177° C (350° F) 





Material properties 


Young's modulus 

j 

El 

. 

GPa 

psi 

155 

22.5 X 10® 

1 

155 

22.5 X 106 

E2 

GPa 

psi 

8.83 

1.28 X 1q6 

5.58 

0.81 X 106 

Poisson's ratio, major 

Vi2 


0.30 

0.23 

Shear modulus 

^12 

GPa 

i 5.10 

0.551 



psi 

o 

X 

• 

o 

0.08 X 106 

Coefficients of thermal 


per °C 

-0.234 X 10~6 

-0.126 X 10-6 

expansion 


per ^ 

-0.13 X 10"® 

-0.07 X 10-6 


0^2 

per °C 

34.0 X 10-6 




per °F 

18.9 X 10-6 










TABLE 2.- Continued 


(a) Concluded 


Proper ty 


Value at - 


Room temperature 177^ C (350^ F) 


Material allowables 


Strength component in fiber 
direction, tensile and 

xt 

GPa 

psi 

1.11 

0.161 X 106 

1 .03 

0.150 X 10® 


cr o o X V w 

Xc 

MPa 

psi 

-972 

-0.141 X 10^ 

" 

1 -848 

1 -0.123 X 10^ 



Strength component in direction MPa 

normal to fibers, tensile and psi 

compressive 

Yc MPa 


35.8 

5.19 X 1q3 
-170 

-24.7 X 1q3 


21 .0 

3.04 X 103 
-119 

-17.3 X 10^ 


Shear strength 


57.9 

8.40 X io3 


25.9 

3.76 X 103 


Applied loads 


-700 

-4000 


ly kN/m 

Ibf/in. 


xy kN/m 

Ibf/ in. 




















TABLE 2.- Concluded 


(b) Thermal properties 


Property 

Symbol 

Units 

Material 

Graphite/Epoxy 

Insulation 

Mass density 

p 

kg/m^ 

1550 

144 



lb/in3 

56 X 10"3 

5.2 X 10”3 

Thermal conductivity 

k 

W/m-°C 

1 .64 

0.18 



Btu-in/sec-ft2-°F 

3.16 X 10‘3 

0.35 X 10"3 

Heat capacity 

c 

J/kg-°C 

962 




Btu/lb-°F 

0.23 







TABLE 3.- FINAL DESIGNS FOR TITANIUM PLATE WITH ALUMINUM BARS 


Prcx:edure 

Final 

mass 

Design variables 

Maximum 

temperature 

error 

CPU^ per 
iteration, 
sec 

CPU^ for total 
optimization, 
sec 

kg 

lb 

Symbol 

SI 

Units 

U.S. Customary 
Units 

OC 

Op 

Full analysis 

5.35 

11.79 

ti 

21 . 23 mm 

0.836 in. 

0 

0 

1 .99 

167 





6.91 

.272 








^3 

4.70 

.185 








t4 

2.62 

.103 








ai 

901 .9 mm^ 

1.398 in2 








32 

2491 .0 

3.861 








33 

314.8 

.488 



1 


Reciprocal expansion 

1 

1 6.41 

1 

1 14.14 

i ’'1 

28 . 8 5 mm 

, 1 .136 in. 

49.9 

1 89.8 

1 

1 1.12 

100 


1 


^2 

8.70 

I .342 

i 


1 





i t3 

4.66 

: .183 

\ 

i 

' 





t4 

2.76 

.109 

‘ 


i 





i ai 

1036.8 mm^ 

1.607 in^ 








32 

; 2507.7 

1 3.887 







; 

1 

j 342.6 

.531 

1 






^Central processing unit time on CDC Cyber 173; FTN compiler; OPT = 2 







TABLE 4.- FINAL DESIGNS OF INSULATED COMPOSITE PANEL 


[initial design; t-| 
^45 ~ 1.52 mm. 


= 101.6 mm, to = 0.76 nm, 
and tgo =0.76 mm] 


Procedure 

Final 

mass, 

kg/m^ 

(lb/ft2) 

Design variables, mm (in.) 

Maximum 

temperature 

CPU^ per 
analysis, 
sec 

Total 

CPU,^ 

sec 

tl 

to 

t45 

tgo 

error , 
°C (°F) 

Full analysis 

19.94 

(4.084) 

121 .4 
(4.781) 

0.848 

(0.0334) 

0.175 

(0.0069) 

0.546 

(0.0215) 

0 

(0} 

11.00 

230.0 

Direct expansion 

' 

19.98 

(4.093) 

121 .7 
(4.792) 

0.859 

(0.0338) 

0.173 

(0.0068) 

0.541 

(0.0213) 

9.17 

(16.5) 

0.49 

59.1 

Reciprocal expansion 

20.70 

(4.239) 

113.5 

(4.468) 

0.724 

(0.0285) 

1.384 

(0.0545) 

0.688 

(0.0271) 

3.67 

(6.6) 

0.49 

63.0 

Direct expansion 
plus critical 
time points 

19.97 

(4.090) 

121 .7 
(4.790) 

0.853 

(0.0336) 

0.163 

(0.0064) 

0.528 

(0.0208) 

9.50 

(17.1) 

0.23 

28.9 

Reciprocal expansion 
plus critical time 
poi nts 

20.68 

(4.235) 

113.6 

(4.472) 

0.673 

(0.0265) 

1.422 

(0.056) 

0.681 

(0.0268) 

3.72 

(6.7) 

0.23 

28.7 


^Central processing unit time on CDC Cyber 173; FIN compiler, OPT = 2. 




TABLE 5.- FINAL DESIGNS OF INSULATED COMPOSITE PANEL 


jlnitial design: ti = 203.2 iraii and tg = t 45 = tgg = 0.51 


Procedure 

Final 

mass, 

kg/m^ 

(lb/ft2) 

Design variables, mm (in.) 

Maximum 

temperature 

error, 

°C (°F) 

CPu'' per 
analysis, 
sec 

Total 

CPU,'' 

sec 

tl 

to 

t45 

t90 

Full analysis 

19.93 

(4.081) 

121 .5 
(4.783) 

0.848 

(0.0334) 

0.163 

(0.0064) 

0.549 

(0.0216) 

0 

(0) 

11.00 

208.0 

Direct expansion 

13.91 

(2.849) 

80.6 

(3.175) 

0.864 

(0.0340) 

0.102 

(0.0040) 

0.508 

(0.0200) 

292 

2(525) 

0.49 

53.6 

Reciprocal expansion 

20.73 

(4.246) 

127.6 

(5.023) 

0.833 

(0.0328) 

0.124 

(0.0049) 

0.551 

(0.0217) 

66.7 

! (120) 



0.49 58.6 

Direct expansion 
plus critical 
points 

14.49 

(2.968) 

79.9 

(3.145) 

0.777 

(0.0306) 

0.480 

(0.0189) 

0.663 275 0.23 21.5 

(0.0261) J 2(495) 

1 

Reciprocal expansion 
plus critical 
points 

20.67 
1 (4.234) 

127.4 

(5.014) 

0.836 

(0.0329) 

0.114 

(0.0045) 

0.538 ! 68.3 0.23 21.4 

•(0.0212) (1 23) ■ 


^Central processing unit time on CDC Cyber 173; FTN compiler, OPT = 2. 
^Infeasible design. 






T = -18°( 


161 

121 


Percent change 


in T, 


max 




design 

Final 

design 

Percen- 

. S. Cust. 

S.I. unit 

U.S. Cust. 

different 

.516 in. 

21.03 mm 

0.828 in. 

+60.5 

.460 

6.92 

.272 

1 

o 

00 

.164 

4.69 

.185 

+12.2 

.128 

2.62 

"".103 

-19.5 

. 008 in 2 

901.3 mm2 

1.397 in2 

+38.6 

TmY 

2489.7 

3.859 

+ 5.9 

00 

316.1 

.490 

-24.4 


/ 

/ 

/ 

/ 

/ 


Full analysis 
Direct expansion 
Reciprocal expansion 










Figure 6.- Quality of approximations in estimating temperature dependence on insulation 

thickness; T = 2000 sec. 
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Insulation thickness, t, mm 


Figure Variation with insulation thickness of maximum laminate temperature 
and laminate temperature at T = 3200 sec for insulated composite panel. 
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Iteration 

Figure TO,- Performance of Taylor series approximations in optimization 
of titanium plate with aluminum bars. 
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